THE GAMBIER MAPPING 



B. Grammaticos 
LPN, Universite Paris VII 
Tour 24-14, B^'^^etage 
75251 Paris, France 

A. Ramani 
CPT, Ecole Polytechnique 
CNRS, UPR 14 
91128 Palaiseau, France 



Abstract 

We propose a discrete form for an equation due to Gambier and which belongs to the 
class of the fifty second order equations that possess the Painleve property. In the 
continuous case, the solutions of the Gambier equation is obtained through a system of 
Riccati equations. The same holds true in the discrete case also. We use the singularity 
confinement criterion in order to study the integrability of this new mapping. 



1. Introduction. 

One century ago, Painleve classified all second order of the form w" = f{w', w,z), where 
/ is rational in w' , algebraic in w and analytic in z, characterized by the absence of 
movable singularities [1]. This property, that later came to be known as the Painleve 
property, is considered today as an indication of integrability [2]. Integrability was, of 
course, what Painleve had in mind when he undertook his classification. The analysis 
presented by Painleve, however, was not complete. Although he obtained most inter- 
esting results, discovering the equations that define the "Painleve transcendents", his 
analysis was plagued by an oversight. In particular he missed a whole class of equations 
the dominant part of which reads: 

,, n - 1 x'^ 
x" = + . . . 

n X 

As consequence he did not obtain the general forms of the transcendents Piv, Pv and 
Pvi- The gap in the Painleve classification was filled shortly afterwards by Gambler 
who analyzed the problem carefully [3] and presented the complete list of all 50 second 
order equations that possess the Painleve property. 

Foremost among these new equations obtained by Gambler is equation XXVII in 
his classification: 

n — 1 x'^ , , , n — 2 x' no?' o nia' — ah) o „ , 1 , , 

x" = \-axx'^-hx' Ht -^x'^^— 7r-^x^ + fx-h (1.1) 

n X n X (n + 2)^ n+2 nx 

This equation is probably the most complicated that Gambler had to study. It is not 

one of the transcendental equations, neither is it integrable by quadratures. Rather, 

it belongs to the class of equations that can be integrated through linearization. The 

aim of the present paper is to study the discrete analog of this equation of Gambler. In 

previous publications [4,5], we have presented detailed studies of the discrete Painleve 

equations, while the case of linearizable mappings was studied in [6]. However, in the 

latter, we have limited ourselves to mappings that are linearizable through the discrete 

equivalent of the Cole-Hopf transformation. The Gambler equation needs a different 

treatment. It is in fact represented by two Riccati equations 'in cascade'. In what 

follows, we shall present a brief review of the continuous case. This is motivated by 

the fact that Ince's book [7] (the standard reference on Painleve equations) does not 

treat the Gambler equation in any detail. No indication as to how this equation can be 

integrated is given, and its particular cases are not even hinted at. Once the continuous 

equation is properly analyzed, its discretization follows in a most logical way. 



2. A BRIEF RECALL OF THE CONTINUOUS GAMBIER EQUATION. 

The very essence of the Gambier equation is that it is obtained as a cascade of two 
Riccati equations [3] . One starts from a first Riccati: 

y' = -y'' + hy + c (2.1) 

where 6, c are functions of the independent variable z, and couple y to x through a 
second Riccati: 

x' = ax"^ + nyx + a (2-2) 

Here, a is a function of 2;, cr a constant that can be scaled to 1 unless it happens to 
be 0, and n is an integer. (As we shall see in what follows, this last point is a first 
requirement for the absence of critical singularities). The quantities a and a in (2.2) 
are related through a duality. Indeed, replacing x by —1/x we find that the equation 
retains the same form with o and cr interchanged (and n —n). Then if o 7^ 0, the new 
a can be set to 1 (through x ^ x/a, and up to a translation of y: the new a becomes 
aa instead of just a). In order to study the movable singularities of the coupled Riccati 
system we start from the observation that from (2.1) the dominant behaviour of y can 
only be y ~ 1/(2 — zq). The next terms in the expansion of y can be easily obtained, 
and involve the functions b, c and their derivatives. In order to study the structure of 
the singularities of (2.2), we first remark that since the latter is a Riccati, its movable 
singularities are poles. However, (2.2) also has singularities that are due to the singular 
behaviour of the coefficients of the r.h.s of (2.2), namely y. Now, the locations of the 
singularities of the coefficients are 'fixed' as far as (2.2) is concerned. However, from 
the point of view of the full system(2.1-2), these singularities are movable and thus 
should be studied. The 'fixed' character reflects itself in the fact -1 is not a resonance. 
(The terms 'resonance' is used here following the ARS terminology [8] and means the 
order, in the expansion, where a free coefficient enters. A resonance -1 is related to the 
arbitrariness of the location of the singularity, and is thus absent when the location of 
the singularity is determined from the 'outside' rather than by the initial conditions). 
Because of the pole in y, x has a singular expansion with a resonance different from -1 
which may introduce a compatibility condition to be satisfied. 

Before proceeding to the examination of the general case (2.2) let us briefly consider 
the case o = 0, whereupon (2.2) becomes linear. If, moreover, we take cr=0, we flnd 
that the behaviour of x is obtained through x' /x = ny and thus x ~ A{z — zq)"'- In 
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this case, since n is integer, (2.2) always has the Painleve property and the functions b 
and c are free. Next, we turn to the case of the general linear (2.2): x' = nyx + a with 
cr 7^ 0. Let us examine a singularity of the form 

X ~ X{z — zq) (2-3) 

We find X = a/{l — n) unless n = 1. (This last case is excluded as being of non-Painleve 
type. Indeed, one easily finds that for cr = 1, the singular expansion has a logarithmic 
branching point of the form x (z — zq) log(z — ^o)-) Thus, excluding the case n = 1, 
we study just the singular behaviour (2.3). The resonance associated to the singularity 
(2.3) turns out to be exactly n — 1, which explains the requirement that n be an integer. 
For every value of n > 1 we expect a resonance condition (with complexity increasing 
with n). For any n, a sufficient condition is cr = (which, as we have seen above, is also 
necessary for n = 1) and we have a further possibility which, for the first few orders of 
n, writes: 

n = 2 b = 
n = 3 26^ - 6' - c = (2.4) 

n = 4 66^ - 7bb' + b" - 86c + 2c' = 

For n < there is no further constraint. In fact the sufficient condition is a = which 
we have assumed in this paragraph. 

We now turn to the case of the full Riccati (2.3) with a ^ 0. In this case, thanks 
to the duality, we can assume cr = 1 (otherwise interchanging a and cr we are back to 
the previous case). Again, only the singularity due to y can lead to trouble. Rewriting 
(2.2) as x' /x = ax + ny + a/x for y = l/{z — zq) + ... we remark that unless n = ±1 a 
behaviour of the form x ~ (z — zq)"' is impossible when aa ^ 0. For n = 1, a logarithmic 
leading behaviour will be present for cr 7^ 0, for nonzero a exactly as in the case a = 0. 
Conversely the condition cr = is sufficient for the absence of a critical singularity for 
n = 1 even for nonvanishing a. Similarly,in a dual way, for n = —1 the necessary and 
sufficient condition for the absence of a critical singularity is a = 0, irrespective of the 
value of cr. 

Next we assume n 7^ ±1, in which case it suffices to study the singularities x ^ 

\{z — Zq), (A — a/{l — n)) and x ^ iJ,/{z — zq), (// = — (n + l)/a). The first singular 
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behaviour (x ~ X{z — zq)) has a resonance at n — 1, which is negative for n < 1 and 
thus does not introduce any further condition. For n > 1, the resonance condition can 
be studied at least for the first few values of n. We find that a = suffices for the 
resonance condition to be satisfied even for a ^ 0. However, if we demand u 0, we 
find the further possibilities: 

n = 2 6 = 

n = 3 26^ - 6' - c + aa/2 = (2.5) 

n = 4 186^ - 2166' + 36" - 246c + 6c' + Saba - 2a'a = 

The second singular behaviour, x ^ iJ,/{z — zq), has a resonance at — 1 — n. Thus 
for n > this resonance is negative and does not introduce any further condition, while 
for n < a compatibility condition must be satisfied. We find that for every case n < 0, 
a = is a sufficient condition for the absence of critical singularity. (This is not in 
the least astonishing given the duality of a and a). On the other hand if we demand 
a 7^ then a different resonance condition is obtained, at each value of n. For n = — 1, 
whenever o 7^ 0, a logarithmic singularity of the form {z — ZQ)~^log{z — zq) appears 
irrespective of the value of a. For n < —1, we find: 

n = —2 a' — ab = 

n=-3 4a6^ - 2a6' - 2ac + a^a - 6a'6 + 2a" = (2.6) 

n =-4 a(186^-2166'+36"-246c+6c'+8a6(7)+a'(126'+12c-336^-6ac7)+186a" !-3a'"= 

A remark is in order at this point. The analysis presented above was based ex- 

plicitely on the assumption that (2.1) is a Riccati. However, if we take n ^ 00 in 

(2.2) and rescale y (and c) in an appropriate way we find that (2.1) becomes linear. 

y' = by + c. Thus, while on the one hand n ^ 00 in (2.5) leads to a resonance condition 

that should be implemented after an infinite number of steps, on the other hand we do 

not have to consider it. Equation (2.1) being now linear, y does not have any singularity 

that could interact with (2.2). The equation for x reads in that case: 

x'^ / (J \ 

x" = h ( ax + 6 ) x' + (a' — ab^x^ + cx — ab (2-7) 

X \ xf 

which is a form of the equation XIV in the Painleve-Gambier classification [3,7] (not 
the canonical one but an equivalent, generic one). We recall that the canonical form of 
the latter is: 

11 

!(;"=— + (qw + —)w' + ((nP' - r' (2.8) 
w w 
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Equation (2.7) can be integrated to a Riccati: 



x' = ax^ + qx + (j (2.9) 
where g is given by {qcj))' = c(f) and is an auxiliary function related to b through 

b = -070. 

3. The discrete analog of the Gambier equation. 

In order to derive the discrete analog of the Gambier equation we start with the obser- 
vation that in the continuous case we have two Riccati's in cascade. It is well known [6] 
that the discrete equivalent of the Riccati equation is a homographic mapping. Guided 
by this fact we can propose the following form for the discrete Gambier system: 

-by + c 

» = iTT 

coupled to 

X = (3.2) 

1 — ax 

where x = Xm+i, y = Vm+i, and the x, y are understood as Xm, ym- Appropriate 
scalings have been made in the two equations in order to bring them to the forms (3.1- 
2). All four a, 6, c, d are functions of the discrete variable m while cr is a constant that 
can be set to 1 unless it is zero. Let us remark that d may not vanish identically, since 
in that case the mapping (3.2) would decouple. (From (3.1-2), it is straightforward to 
obtain a 3-point mapping for x alone, but the analysis is simpler if we deal with both y 
and x). 

The main tool for the investigation of the integrability of discrete systems is the 
singularity confinement criterion we introduced in [9]. It is based on a conjecture that 
the movable (i.e. initial condition dependent) singularities of integrable mappings dis- 
appear after some iteration steps. It will be used here to investigate the integrability of 
the coupled mapping (3.1-2). 

A first remark before implementing the singularity confinement algorithm is that 
the singularities of a Riccati mapping are automatically confined. Indeed, if we start 
from w = (aw + P)/{'yw + S) and assume that at some step w = —S/^, we find that 
w diverges but w and all subsequent tu's are finite. Thus, the intrinsic singularities of 
(3.2) do not play any role. However, the singularities due to y (obtained from (3.1)) 
may cause problems at the level of (3.2). Two types of singularities may appear. Either 
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y diverges (related to the singularity pattern for y: {. . . , —1, oo, 6, . . .}) or y takes the 
value y = —aa/d. In the latter case, we find x = a irrespective of x, and thus x loses 
one degree of freedom. On the other hand, once we enter a singularity there is no 
way for x to leave it unless y assumes again a singular value (after a certain number 
of steps). Thus the first requirement for confinement would be for y to become again 
either oo or equal to —aa/d. However, if y were to take the value cxo again some steps 
after first taking it, it would take it periodically and the singularity would then also be 
periodic. This is contrary to the requirement that the singularity be movable: a periodic 
singularity (with fixed period) is 'fixed' in our terminology. The same applies to the 
singularity y = —aa/d. Thus the only singularity pattern for y that are acceptable 
are: either start with —aa/d and reach oo after N steps, or start with oo and reach 
—aa/d after N' steps. (Clearly, these two singularity patterns are mutually exclusive: 
one cannot combine them since in that case oo would be followed by oo after N + N' 
steps and thus the singularity would again be periodic). So the first condition is that 
y take the value —aa/d N steps before (or A^' steps after) taking the value oo. The 
equivalent of this requirement in the continuous case is that the resonance be integer. 
We see here that this condition becomes here more complicated. In fact, the first few 
read: 

N=l -aa/d=-l 

N = 2 -aa/d= -P^ (3.3) 
' b+1 ^ ^ 

N = 3 —aa/d = — = 

{b+l)b + c+l 

In what follows, we shall label these cases ior N > while the cases where the value 
—aa/d appears N' steps after oo, such as: 

N' = 1 - aa/d = b (3.4) 

bb + c 

N' = 2 - aa d = 

wiU be labelled by AT = -N' < 0. 

While one can compute easily more such conditions, it is not possible to give a 
general expression, like the "n must be integer" of the continuous case. 

Once y hits a second singularity there is a possibility for x to recover its lost degree 
of freedom through some indeterminate form (like O.oo or 0/0). Thus when y becomes 



infinite {N steps after assuming the value —aa/d) we must simultaneously have a; = 0. 
Similarly when y assumes the value —aa/d {N' steps after being infinite) x must be 
equal to 1/a so as to lead to 0/0 for x. We can study the first few cases and obtain 
the condition for the confinement of the singularity. Let us start with N < {). For 
N — —1, we enter the singularity through y = oo, and thus x = oo, y = b. We ask 
next that —aa/d = y = b and moreover that this lead to an indeterminate form in the 
calculation of x. The condition in that case is just a = (and this implies b = also, 
as d identically zero is not allowed). 

For N = —2 we obtain x = —db/a and ask that 

= bb + c na 
^+1 d 

The condition for x to result from an indeterminate form is 1 — ax = i.e. adb + a = 0. 
This last condition is automatically satisfied when a = 0. 

In fact, it turns out that for every A'"<0, a = Oisa sufficient condition for con- 
finement. In this case, the second equation of the mapping becomes linear. Obtaining 
the conditions for the next higher N is straightforward. For N = —3, for instance, we 

find: _ _ _ 

= bbb + bc + bc + c aa 

^ bb + b + c + 1 I 

leading to 

{b+l){a^a -a-ndb) -^ldb{bb + c) = (3.6) 

The second possibility is to enter the singularity through y = —aa/d, i.e. the N > 
cases. For N = 1, we enter the singularity through y = —aa/d which gives x = a and 
we exit through y = oo, which implies aa/d = 1. In order to have an indeterminate 
form at the level of x we must have xy = O.oo and thus cr = a condition clearly 
incompatible with aa/d = 1. Thus (3.1-2) as written cannot confine for A?" = 1 . In 
fact, the confinement condition would be o" = provided the denominator of (3.1) is 
just y instead of y + 1 (which is dual to = — 1 with a = b = 0). However, for > 1 
we will find that it is possible to recover the lost degree of freedom. For = 2, we ask 
that x be zero and it turns out that this is possible if either a — or d = 1. In fact, the 
condition cr = is sufficient for all cases A'^ > 1 (provided, of course, that the condition 
for AT, i.e. "y becomes infinite after N steps", holds) For N = 3, we write the condition 
(based on the assumption a ^ 0): 

{b+l){aa + d- 1) - {l + c)dd = (3.7) 
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Conditions for the next higher A'"'s can easily be obtained. 

A question that naturally comes to the mind is whether the continuous limit of 
system (3.1-2) is indeed the coupled-Riccati original system of Gambler. The answer 
appears trivially affirmative at first sight: in fact, system (3.1-2) was constructed as 
a direct discretization of Gambler's system. However, a subtle question remains. In 
equation (2.2) of Gambler, n appears explicitely: what is its equivalent in the discrete 
case? On the one hand, it is always possible to present a formal limit of (3.1-2) where n 
is related to d; in fact we have d = n — 1. However, in the previous paragraph we have 
argued that what plays the role of the position of the resonance in the discrete case is 
a complicated relation. This condition alone is not sufficient in order to ensure a form 
like (2.2) at the continuous limit with n appearing naturally. (In particular, it does not 
prevent d from being non constant, which at the continuous limit would have given a 
n function of z). Only when one takes into account both confinement conditions does 
system (3.1-2) turn out to be of the precise form (2.1-2) at the continuous limit for the 
corresponding n. Moreover the 'continuous' n is precisely the 'discrete' introduced by 
the confinement. (Note, for instance, that the second confinement condition for N — 2 
is, precisely, d—1 i.e. n = 2). 

The discussion presented above was based on the form (2.1) for the first Riccati. 
However one should also consider another case that cannot be obtained as a limit of 
(2.1) with the scaling we chose, namely: 



Once the denominator has been simplified a new scaling is possible and we rewrite (3.8) 
as: 



The analysis of the coupled system (3.9) and (2.2) proceeds along the same lines as 
before. One obtains simpler discrete forms of the Gambler equation with only one free 
function instead of two in the generic case. As we have already remarked, this is the 
only case that confines for A?" = 1, the condition being still cr = 0. 
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y = b+- 

y 



(3.8) 



y = b+- 

y 



(3.9) 
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4. Discrete systems related to the Gambier mapping. 

In the previous section, we have examined the mappings consisting in two Riccati's 
in cascade. Moreover, the particular case where the second Riccati becomes a hnear 
equation was included in our study through a = and, because of the duality of a and 
cr, through cr = as well. Finally the case where the first Riccati reduces toy = b + c/y 
was briefly considered in section 3 since it does not alter qualitatively the conclusions 
reached with the full (3.1). In this section, we are going to consider two particular cases 
of the Gambier mapping. The first case corresponds to equation (3.1) becoming linear. 
We thus have the mapping: 

y^y + c (4.1) 

X = (4.2) 

1 — ax 

(where by the appropriate scaling we have taken 6=1). 

Clearly, this mapping is a discrete form of the equation P14 encountered in section 
2. This can be readily seen when we consider its continuous limit. Putting x = w/e, 
d = eq{z), a = e'^p{z), c = —s + er{z)/d where s = q' /q, we obtain indeed, at the limit 

w" = h (pw + s ] w' + {p' — sp)w^ + (r + s^)qw — ah (4-3) 

w \ w/ 

i.e. exactly (2.7), i.e. P14, in noncanonical form. 

We thus expect system (4.1-2) to be integrable without constraint. What does 

singularity confinement tell us on this point? Before proceeding to the actual calculation, 

let us try to find the answer by analogy. The case of a linear -|-Riccati system was 

obtained in the continuous case by the limit n — > cxo. Since the number N of steps 

before confinement is related to the position n of the resonance in the continuous case 

(in fact n = N) we expect confinement to be obtained only after an infinite number 

of steps! This turns out indeed to be the case. Given the structure of the mapping 

(4.1-2) the only singularity can he y = —aa/d (where (4.2) loses one degree of freedom). 

However, the lost degree of freedom cannot be recovered in a finite number of steps 

since we need y = 00. Thus the mapping has a very peculiar behaviour: it can only 

confine after an infinite number of iterations steps. More on this subtle point will be 

given in the next section. 
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Let us now turn to a second particular case of the Gambier mapping. Suppose 
now that the second equation (3.2) of the mapping has a denominator equal to —ax 
instead oil — ax. To avoid a trivial reduction we must here take cr 7^ 0. Then (after a 
rescaling and a redefinition of the parameters) the system becomes: 

V = ^ (4.4) 
x=^^M^ (4.5) 

X 

Here also we have a peculiar behaviour with respect to singularity confinement. Clearly, 
(4.5) loses one degree of freedom when y = 00. Moreover, there is no possibility to 
recover this lost degree of freedom for any value of y. Thus the singularity never confines. 
As we shall argue in the next section this is a genuinely nonconfined singularity. Our 
argument can be further strengthened by considering the continuous limit of (4.4-5). 
We put a: = l + ew,b= -1 - e - 6^(5 - 1/2)), c = b + e^{-l/4 + e{2q - l)/4 + e'^p) and 
take just d = —2 + e, the latter being clearly a restriction which is sufficient in order to 
prove our point. We find 

w" = -2w'w - (w' + w'^f - 2q{w' + w'^)-q + 4p+^ (4.6) 

This last equation is manifestly of non-Painleve type (because of the {w' + w'^)'^ term). 
Although (4.4-5) is not confining in general, there exists a particular case where one 
can recover confinement, namely when (4.4) becomes linear: 

y = y + c (4.7) 

^=^^ (4.8) 

X 

In this case, ther is no way for y to become singular and thus the mapping has only 
confined singularities. The continuous limit is compatible with our conclusion. Indeed, 
putting X = 1 + ew, d — d{z), c = —2ed'/d'^ + c'q/d we find at the continuous limit 
e — > the equation: 

d' 

w" = -2w'w -—(«;' + w^) - q (4.9) 

The latter is (a non canonical form of) equation P5 in the Painleve-Gambier list [3,7] 
and is obtained precisely as the derivative of a Riccati. (Note that by P5 we mean here 
the fifth equation in the list and not the fifth Painleve transcendent, which we usually 
denote by Py)- 

11 



Thus, the confining subcases of the Gambier mapping obtained in this section 
are all variants of the 'discrete derivative' (with respect to one of the parameters) 
of a homographic mapping. This means that we can start with a mapping of the 
form w = ^j^, solve for one of the parameters, take the discrete derivative of this 
parameter and assign an m-dependent value to it. Given the duality of the homographic 
mapping under the transformation w 1/w, only two different types of such mappings, 
corresponding to derivatives with respect to either a or P, exist. 

5. A DISCUSSION OF THE CONFINEMENT PROPERTY. 

In the preceding sections, we have repeatedly used the criterion of the confinement of 
singularities in order to characterize the integrable or not character of the Gambier 
mapping under its various forms. Since the behaviours of the systems of section 4 
were rather atypical as far as the confinement is concerned, we feel that some detailed 
discussion is necessary at this point. 

The confinement property was initially introduced (and related to integrability) 
for rational mappings, the structure of which would lead to vanishing denominators 
and the appearence of divergences [9]. However it was rapidly realized that this is 
not the only manifestation of a singular behaviour. Thus, in [6], we introduced the 
notion of singularity so as to include all cases where the system loses one degree of 
freedom. Confinement of the singularity in this case means recovery of the lost degree 
of freedom through the appearence of an indeterminate form (0/0, O.oo, etc.). One 
important condition (in perfect analogy to the continuous case) is that the singularity 
be 'movable' i.e. initial condition dependent. Let us illustrate this last point in the case 
of the Riccati mapping: 

Wm+l = —J (5.1) 

where ct, /3, 7 and 5 depend in general on m. The mapping variable w can very well 
become infinite for some m, but this (movable) singularity is immediately confined: 
infinity does not play a particular role in a homographic mapping. However, if at certain 
iteration the determinant a5—j3^ of the coefficient matrix vanishes the mapping loses one 
degree of freedom and moreover there is no way that it can recover it in its subsequent 
evolution. However this non-confined singularity is perfectly acceptable since it is fixed 
i.e. its location does not depend on the initial conditions. 

Clearly, the Riccati mapping is too simple to exhibit any interesting behaviour as 
far as the singularity confinement is concerned. As soon as we enlarge our scope and 
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consider systems of two Riccati's in cascade (as in the case of the Gambier mapping) 
the behaviour becomes considerably richer. A first point must be made about periodic 
singularities. They are assimillated to fixed ones since they do not have a specific 
location, that could depend on the initial conditions, but rather they extend everywhere. 

The main difficulty of the Gambier mapping comes from the fact that the dis- 
tinction between fixed and movable singularities is somewhat artificial [10]. A movable 
singularity in the first equation of the system enters in the second one through a coeffi- 
cient and should thus be considered as 'fixed', as far as the second equation is concerned. 
(Let us recall again that the 'movable' singularities of a homographic mapping (5.1) do 
not cause any difficulties). Still its location depends on the initial conditions of the 
system, considered globally, and must, thus, be treated as a movable singularity. 

In section 4 we have encountered two pathological cases. In one of them, we have 
argued that confinement was possible (although after an infinite number of steps) and 
thus that the equation was in principle confining. The second equation was rejected as 
non-confining since there was no possibility whatsoever to confine. We are aware that 
the distinction is subtle and lies on grounds that are not very firm. In the cases at 
hand, we were guided by the underlying structure of the system in order to make our 
claims on the confining or not character of the mapping. However, if we are given a 
general 3-point mapping, there is a priori no easy, algorithmic, way to decide whether 
the mapping confines after an infinite number of steps, or not at all. We hope that just 
as in the continuous case, the occurence of a resonance aX n — oo, compatible with the 
Painleve property, is an exceptional occurence, and that the 3-point mappings confining 
after an infinite number of iterations will be restricted to those of the Gambier family. 

A useful guide as to the confining or not character of the mappings (in such tricky 
cases) is the study of all possible continuous limits. Our conjecture is that confining 
mappings have as continuous limits differential equations that possess the Painleve 
property. Thus if for some mapping we are able to find some non-Painleve continuous 
limit, this is a clear indication as to its non-confining character. However, the absence 
of such non-Painleve limits does not allow us to draw clear conclusions: it may only 
mean that non-Painleve continuous limits do not exist (or that we were not able to find 
them!) Still, the continuous limit can be a 'last stand' argument for some mapping that 
is difficult to interpret. 

Finally, let us point out that all our discussion of Gambler's mapping and, in fact, 

13 



of Gambier's equation as well, is a discussion of the Painleve property rather than plain 
integrability. In a system such as (2.1-2) or (3.1-2) of Riccati equations in cascade, 
we can always solve the first equation for y, obtain y{t) or j/^, and inject it into the 
second equation. The latter can always be written as a linear second order difi"erential 
equation for x, or a linear 3-point mapping. So, in principle, the problem can always be 
solved formally. The difficulty comes when one wishes to actually compute x, in terms 
of contour integrals in the continuous case while y has bad analytic properties. The 
Painleve property guarantees that one can perform the integration and obtain indeed 
the solution of x{t) over the complex t-plane. In the discrete case the difficulty arises 
when one tries to compute Xm in terms of matrix products (the elements of which 
contain y^). When has the wrong properties some of these matrices are singular 
in such a way that degrees of freedom are irretrievably lost. Singularity confinement 
means that these lost degrees of freedom are recovered at some later stage. 

6. Conclusion. 

In this paper we have analysed a mapping consisting in two discrete Riccati's in cascade. 
This mapping is the discrete analogue to Gambier's equation (number XXVII of the 
Painleve- Gambler list). The interesting feature of this discrete system is the fact that 
the singularitites can be confined after an arbitrary number of steps. This is quite 
unusual, since, in most cases encountered up to now, confinement was obtained after a 
small number of iteration steps. This new behaviour, encountered for the first time in the 
Gambler mapping has made possible the refinement of the whole notion of singularity 
confinement. 
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